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Abstract It appears that string-M-theory is the only viable candidate for a com- 
plete theory of matter. It must therefore contain both gravity and QCD. 
What is particularly surprising is the recent conjecture that strongly 
coupled QCD matrix elements can be evaluated though a duality with 
weakly coupled gravity. To date there has been no direct verification of 
this conjecture by Maldacena because of the difficulty of direct strong 
coupling calculations in gauge theories. We report here on some progress 
in evaluating a gauge-invariant correlator in the non-perturbative regime 
in two and three dimensions in SYM theories. The calculations are 
made using supersymmetric discrete light-cone quantization (SDLCQ). 
We consider a Maldacena-type conjecture applied to the near horizon 
geometry of a Dl-brane in the supergravity approximation, solve the 
corresponding N = (8, 8) SYM theory in two dimensions, and evaluate 
the correlator of the stress-energy tensor. Our numerical results sup- 
port the Maldacena conjecture and are within 10-15% of the predicted 
results. We also present a calculation of the stress-energy correlator in 
N = 1 SYM theory in 2-1-1 dimensions. While there is no known duality 
relating this theory to supergravity, the theory does have massless BPS 
states, and the correlator gives important information about the BPS 
wave function in the non-perturbative regime. 
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Introduction 

Recently the conjecture that certain field theories admit concrete re- 
alizations as string theories on particular backgrounds has caused a lot 
of excitement. The original Maldaccna conjecture [1] asserts that the 
J\f = 4 supersymmetric Yang-Mills (SYM) theory in 3-1-1 dimensions 
is equivalent to Type IIB string theory on an AdSs x background. 
However, more recently, other string/field theory correspondences have 
been conjectured. Attempts to rigorously test these conjectures have 
met with only limited success, because our understanding of both sides 
of the correspondences is usually insufficient. The main obstacle is that 
at the point of correspondence we want the curvature of space-time to 
be small in order to use the supergravity approximation to string theory. 
This requires a non-perturbative calculation on the field theory side. We 
use the method, SDLCQ, in the corresponding non-perturbative regime. 
SDLCQ, or Supersymmetric Discretized Light-Cone Quantization, is a 
non-perturbative method for solving bound-state problems that has been 
shown to have excellent convergence properties[2]. 

Aside from our numerical solutions, there has been very little work on 
solving SYM theories using methods that might be described as being 
from first principles. While selected properties of these theories have 
been investigated, one needs the complete solution of the theory to cal- 
culate the correlators. By a "complete solution" we mean the spectrum 
and the wave functions of the theory in some well-defined basis. The 
SYM theories that are needed for the correspondence with supergravity 
and string theory have typically a high degree of supersymmetry and 
therefore a large number of fields. The number of fields significantly 
increases the size of the numerical problem. Therefore, when presenting 
the first calculation of correlators in 2-1-1 dimensions, we consider only 
Af = 1 SYM. 

An important step in these considerations is to find an observable 
that can be computed relatively easily on both sides of a string/field 
theory correspondence. It turns out that the correlation function of 
a gauge invariant operator is a well-behaved object in this sense. We 
chose the stress-energy tensor T^'^ as this operator and will construct 
this observable in the supergravity approximation to string theory and 
perform a non-perturbative SDLCQ calculation of this correlator on the 
field theory side. 

1. String Theory 

String theory contains solitons, the so-called D-branes on which modes 
can propagate as well as in the bulk. While in general these modes 
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Figure 1. Phase diagram of two-dimcusioual M = (8, 8) SYM: the theory flows from 
a CFT in the UV to a conformal a-model in the IR. The SUGRA approximation is 
valid in the intermediate range of distances, 1/ gYM\/N2 < x < ^fN'd qym- 

couple, there exists a limit in which the bulk modes decouple from the 
modes on the D-brane; this is typically a low energy limit. In this limit 
the theory on Dj)-brancs, separated by at most sub-stringy distances 
becomes a supersymmetric S\J (N) Yang-Mills theory. As the D-branes 
carry mass and charge, they can excite gravity modes in the bulk, for 
which in supergravity there exist equivalent solutions. One thus has 
a string/field theory correspondence. Naively, one would think that 
supergravity can only describe the large distance behavior of fields, but 
it turns out that one can trust these solutions as long as the curvature 
is small compared to the string scale. In this sense, the large N limit is 
a valid description [3]. 

The most prominent string/field theory correspondence is the so- 
called Maldacena conjecture[l], which assures that the conformal A/" = 4 
SYM in 3+1 dimensions, is equivalent to a type IIB string theory on 
a AdS^ X background. In the more general case of non-conformal 
theories, it turns out that a black p-brane solution, stretching in p -|- 1 
spacetime dimensions, of supergravity will correspond to a supersym- 
metric Yang-Mills theory in p + 1 dimensions [3] . 

One can test these string/field theory correspondences, if one is able 
to construct and evaluate observables on both the string and the field 
theory regimes. Although this at first seems a hard task, because typi- 
cally the small curvature regime of string theory, where the supergravity 
approximation allows quantitative calculations, falls into the strong cou- 
pling regime of the field theory side. We shall see that we can come up 
with scenarios where we can evaluate the field theory observable non- 
perturbatively. 

1.1. Two-dimensional correlation functions from 



It is instructive to take a closer look on the expected properties of 
M = (8, 8) SYM in two dimensions, before we proceed to technical details 
on the string theory side. In the extreme ultra-violet (UV) this theory 
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is conformally free and has a central charge cjjv = ^c- Perturbation 
theory in turn will be valid for small effective couplings g = gYuV^cX, 
where x is a space coordinate. For large distances, in the far infra- 
red (IR), the theory becomes a conformal cr-model with target space 
{R^)'^" / Snc- The central charge is cir = Nc- It is a bit more involved 
to show that here perturbation theory breaks down when X'^\fN'cl 9YM-, 
see e.g. Ref. [3]. 

The intermediate region, 'VjgYM^n^c < x < y/N^/gyM, where no per- 
turbative field theoretical description is possible, is fortunately exactly 
the region which is accessible to string theory; or rather, to the super- 
gravity (SUGRA) approximation to Type IIB string theory on a special 
background. It is that of the near horizon geometry of a Dl-brane in 
the string frame, which has the metric 

ds'^ = a'pYM ( ^dxf, + + C/dQg.p j 




^0 _ '^'^9ym - 

= — ttq — 9YM, [i-) 



where we defined gvM = ^Tr^^^gyMV^- In the description of the com- 
putation of the two-point function we follow Ref. [4] . The correlator has 
been derived in Ref. [5], being itself a generalization of Refs. [6, 7]. 

First, we need to know the action of the diagonal fluctuations around 
this background to the quadratic order. We would like to use the ana- 
logue of Ref. [8] for our background, Eq. (1), which is not (yet) available 
in the literature. However, we can identify some diagonal fluctuating 
degrees of freedom by following the work on black hole absorption cross- 
sections [9, 10]. One can show that the fluctuations parameterized like 



ds^ = (l + /(x°,f/)+5(x°,[/))5oo(dx°)2 

+ (l + 5/(x°, U) + gix", U)) guidx'f 
+ (l + U) + g{x°, U)) guudU^ 
+ + fix', U) - ^-g{x\ 5nad02 

el' = (l + Zf{x\U)-g{x\U))e't>\ (2) 
satisfy the following equations of motion 

f"{U) = -Lf^u) + ^^f{U) (3) 
9'\U) = -yiU) + ^9{U) + '^g{U). 
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Without loss of generality we have assumed here that these fluctuations 
vary only along the direction of the world volume coordinates, and 
behave like a plane wave. One can interpret a Dl-brane as a black hole 
in nine dimensions. The fields f{U) and g{U) are the the minimal set 
of fixed scalars in this black hole geometry. In ten dimensions, how- 
ever, we see that they are really part of the gravitational fluctuation. 
Consequently, we expect that they are associated with the stress-energy 
tensor in the operator field correspondence of Refs. [6, 7]. In the case of 
the correspondence between A/'=4 SYM field theory and string theory 
on an AdS^ x background, the superconformal symmetry allows for 
the identification of operators and fields in short multiplets [11]. In the 
present case of a Dl-brane, we do not have superconformal invariance, 
and this technique is not applicable. Actually, we expect all fields of the 
theory consistent with the symmetry of a given operator to mix. The 
large distance behavior should then be dominated by the contribution 
with the longest range. The field /(A;°,{7) appears to be the one with 
the longest range since it is the lightest field. 
Eq. (3) for f{U) can be solved explicitly 

f{U) = U~'K,,, (l^fc) , (4) 

where K'^i2{x) is a modified Bessel function. If we take f{U) to be 
the analogue of the minimally coupled scalar, we can construct the flux 
factor 

= hm ^V^<7^^e-2(^-^-)9t/log(/(C/))|t;=t;„ 
— u |_ 

up to a numerical coefficient of order one which we have suppressed. 
We see that the leading non-analytic contribution in fe^ is due to the 
term. Fourier transforming the latter yields 

3 

{0{x)Om = (6) 

This is in line with the discussion at the beginning of this section. We 
expect to deviate from the trivial (1/a:^) scaling behavior of the corre- 
lator at xi = '\^/ gYM^/N~c and X2 = VT^d Qym- This yields the phase 
diagram in Fig. 1. It is interesting to note that the entire hierarchy 
is consistent in the sense of Zamolodchikov's c-theorem, which assures 
that the central charges obey c(x) > c(y), whenever x <y [12]. 
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SYM3 1 D2-branes ^ M2-branes 1 AdS4 x S'^ x 

Figure 2. Phase diagram of thrcc-diincrisiorial M' = (8, 8) SYM: the theory flows 
from a perturbative SYM in the UV to a M-theory on AdS4 x S^. The SUGRA approx- 
imation is vaUd in the intermediate range of distances, I/qym^c < x < I/qym^c^^ ■ 



1.2. D2-branes and three-dimensional SYM 

As stated in the introduction, in an analogous way one can show that 
a system of Z'2-branes corresponds in a certain limit to a Yang-Mills 
theory in three dimensions. It is again a A/" = (8, 8) super symmetric 
theory. Unfortunately, an observable like the correlation function of the 
stress-energy tensor has not yet been calculated for this theory. However, 
there are encouraging results both on the string and on the field theory 
side of the correspondence [3, 13]. 

We describe the phase-diagram of three-dimensional SYM with 16 
supercharges here, following Ref. [3]. Later we will present a non- 
perturbative field theory calculation within the SDLCQ framework. The 
latter calculation is, however, of a theory with an A/" = 1 supersymme- 
try. This theory might nevertheless share some features with the full 
M = (8, 8) theory, cf. also the results of two-dimensional SYM with 
different supersymmetries in Sec. 2. 

It can be argued that the theory has to be described by different de- 
grees of freedom at difi'erent energy scales. At large Nc, one can use 
perturbation theory of SYM(2-|-1) in the far ultra-violet, i.e. at small 
distances. The supergravity solution, which in this regime is an ap- 
proximation to type IIA string theory with D2-branes, can be trusted 
at intermediate distances r, l/gyM^c < r < l/gyj^N^^^ . It has been 
conjectured that for large distances, r > l/(7y^A^^/^, an M-theory de- 
scription is appropriate, while in the far infrared, r 3> l/gy^, this theory 
is equivalent to M-theory on an AdS^ x S''' background, dual to a CFT 
with an 6*0(8) R-symmetry. This picture is compiled in Fig. 2. 

2. Field theory correlators and SDLCQ 

Discretized Light-Cone Quantization (DLCQ) preserves supersymme- 
try at every stage of the calculation if the supercharge rather than the 
Hamiltonian is diagonalized [14, 15]. The framework of super symmetric 
DLCQ (SDLCQ) allows one to use the advantages of light-cone quantiza- 
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tion {e.g. a simpler vacuum) together with the excellent renormalization 
properties guaranteed by supersymmetry. 

The technique of (S)DLCQ was reviewed in Ref. [16], so we can be 
brief here. The basic idea of light-cone quantization is to parameterize 
space-time using light-cone coordinates 



and to quantize the theory making play the role of time. In the 
discrete light-cone approach, we require the momentum p- = p'^ along 
the direction to take on discrete values in units of p'^/K where p+ 
is the conserved total momentum of the system. The integer K is the 
so-called harmonic resolution, and plays the role of a discretization pa- 
rameter. One can think of this discretization as a consequence of com- 
pactifying the x" coordinate on a circle with a period 2L = 2itK/p'^ . 
The advantage of discrctizing on the light cone is the fact that the di- 
mension of the Hilbert space becomes finite. Therefore, the Hamiltonian 
is a finite-dimensional matrix, and its dynamics can be solved explicitly. 
In SDLCQ one makes the DLCQ approximation to the supercharges Q*. 
Surprisingly, also the discrete representations of Qi satisfy the supersym- 
metry algebra. Therefore SDLCQ enjoys the improved renormalization 
properties of supersymmetric theories. To recover the continuum result, 
K has to go to infinity. We finds is that SDLCQ usually converges much 
faster than the naive DLCQ. 

In the three-dimensional case we also discretize the transverse momen- 
tum along the direction x"*-; however, it is treated in a fundamentally 
different way. The transverse resolution is T, and we think of the theory 
as being compactified on a transverse circle of length I. Therefore, the 
transverse momentum is cut off at rk.2iTT/l and discretized in units of 
2tt/1. Removal of this transverse momentum cutoff therefore corresponds 
to taking the transverse resolution T to infinity. 

Let us now review these ideas in the context of a specific super-Yang- 
Mills (SYM) theory. Actually, it turns out that the two-dimensional 
SYM is essentially 'included' in the three-dimensional case [17], in the 
sense that in the weak coupling limit the spectrum of the three-dimension- 
al theory is that of the lower dimensional theory. We therefore describe 
here only the 'more general' three-dimensional theory and hint at the 
differences and changes that are to make to recover the two-dimensional 
theory. We start with 2-1-1 dimensional N = 1 super- Yang- Mills theory 
[18] defined on a space-time with one transverse dimension compactified 




(7) 



on a circle. The action is 



(8) 



After introducing the Ught-cone coordinates = -^{x'^ ^x^), decom- 
posing the spinor * in terms of chiral projections 



1 + 7" 
2V4 



X 



1-7^ 



(9) 



and choosing the Ught-cone gauge = 0, we obtain the action in the 
form 



= J dx'^dx j 



dx I tr 



(10) 



A simpUfication of the hght-cone gauge is that the non-dynamical fields 
A~ and x "^^Y be explicitly solved from their Euler-Lagrange equations 
of motion 



A- 



9YM 



J 



^^^(i[</-,a_0] + 2V^V), x = -^^±V'- (11) 



These expressions may be used to express any operator in terms of 
the physical degrees of freedom only. In particular, the light-cone en- 
ergy, P~ , and momentum operators, , corresponding to transla- 
tion invariance in each of the coordinates x^ and x_\_ may be calculated 
explicitly as 



P+ = 
P- = 

P± = 



dx I tr 



dx 



dx 



dx±ti 



{d-(^f + 111)3-11) 



(12) 
(13) 
(14) 



The light-cone supercharge in this theory is a two-component Majorana 
spinor, and may be conveniently decomposed in terms of its chiral pro- 
jections 



Q+ = 2^/^ J dx- j^dx±tT[(f)d-i/)-il)d-(l)], 



(15) 



Q' 



23/4 



dx I tr 



2d±(f)il; + gvM (i[0, 5-0] + 2V'V) 
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The action (9) gives the following canonical (anti-)commutation relations 
for propagating fields for large Nc at equal x"*": 

[(f)ij{x~,x±),d-(l)ki{y~,y±)\ = {'>Pij{x' ,x±),'il^ki{y'' ,y±)} (16) 

= 2'^'^^'" ~ y~)S{x± - y±)Sii5jk. 

Using these relations one can check the supcrsymmetry algebra 

{Q^,Q^+} = 2V2P^, {Q+,Q-} = -iP±. (17) 



In solving for mass eigenstates, we will consider only states which 
have vanishing transverse momentum, which is possible since the total 
transverse momentum operator is kincmatical. Strictly speaking, on a 
transverse cylinder, there are separate sectors with total transverse mo- 
menta 2TrN±/L; we consider only one of them, N± = 0. On such states, 
the light-cone supercharges Q"*" and Q~ anti-commutc with each other, 
and the supcrsymmetry algebra is equivalent to the J\f = (1, 1) super- 
symmetry of the dimensionally reduced (i.e., two-dimensional) theory 
[14]. Moreover, in the Pj_ = sector, the mass squared operator is 
given by = 2P+p-. 

As we mentioned earlier, in order to render the bound-state equations 
numerically tractable, the transverse momenta of partons must be trun- 
cated. First, we introduce the Fourier expansion for the fields (/) and ip, 
where the transverse space-time coordinate x± is periodically identified 



(l)ij{0,x ,x_i_) = 



1 



E 



ipij{0,x ,x±) 



bij {k~^ , 



n 



1 „—ik'^x 



E 



V2k+ 
dk+ 



(18) 



(19) 



''^+bUk+,n^)e 



ik+x-+i'- 



Substituting these into the (anti-)commutators (16), one finds 
aij{p'^,n±),al^{q'^,m±) = [bij{p'^ ,n±),b],^{q+ ,m±)^ 

The supercharges then take the following form: 



Q' 



i2V4 J2 



dkVk 



(20) 



10 

X ^blj{k,n'^)aij{k,n'^) — a\j{k,n'^)bij{k,n^ 



X ^a\j{k,n-^)hij{k,n-^) — h\j{k,n-^)aij{k,n-^)\ + 
gym I Q;/,^Q;^2cifc3<5(/ci + fcs - A:3)<^„x+„x,„x | 

bljih, ni)aik{ki,ni)akjik2,n2 



-alk{ki,ni)blj{k2,n^)aij{k3,n^)] 
^i)aij{k2,ni)aij{k3,ni) 

-aUk^, n^)bik{ki)akj{k2, n^)] 



1^kxk2 


A;3 


1 


fci + ^3 




k2 


1 


k2 + A;3 







+ ^ - •^)[4(^i'«^)4i(fc2,n^)5i,(fc3,n^ 



+6|j(A;3,n3 )6ife(A;i,n5^)6fej(A;2,n2 )] | • 



We now perform the truncation procedure; namely, in all sums over 
the transverse momenta nf appearing in the above expressions for the 
supercharges, we restrict summation to the following allowed momentum 
modes: nf = 0, ±1... ± T. Note that this prescription is symmetric, in 
the sense that there are as many positive modes as there are negative 
ones. In this way we retain parity symmetry in the transverse direction. 
The longitudinal momenta ki = n^Tr/L are restricted by the longitudinal 
resolution according to K = ^^ni. 

The two-dimensional supercharges are essentially recovered, when we 
put n±_ to zero. In particular, the first term of the supercharge Q~ , 
Eq. (21), is absent in this case. Additionally, we have to adjust the 
normalization constants in front of the expressions for the supercharges. 

2.1. Two dimensional correlators 

Using SDLCQ, we can reproduce the SUGRA scaling relation, Eq. (6), 
fix the numerical coefficient, and calculate the cross-over behavior at 
l/^YM-v/^ < X < y/Nc/gYM- To exclude subtleties, nota bene issues of 
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zero modes, we checked our results against the free fermion and the 't 
Hooft model and found consistent results. 

Let us now focus on the theory in two dimensions. We would like to 
compute a general expression for the correlator of the form F{x~ ,x^) = 
{0{x~ , x'^)O{0,0)) . In DLCQ one fixes the total momentum in the x~ 
direction, and it is natural to compute the Fourier transform and express 
it in a spectrally decomposed form 



F(P_,x+) 



2L 



(O(P_,x+)O(-P_,0)) 



Y-{0\O{P.)\n)e-^n-Un\Oi-P-,0m ■ (22) 



The form of the correlation function in position space is then recovered 

by Fourier transforming with respect to P = Kn/L. We can continue 
to Euclidean space by taking r = V2x+x~ 
correlator of the stress-energy tensor is 



L 



to be real. The result for the 



F{x-,x+) = Yl -{n\T++{-K)\0) 



Mi 



K4 (m„V2 




(23) 



where is a mass eigenvalue and K4{x) is the modified Bessel function 
of order 4. Note that this quantity depends on the harmonic resolu- 
tion K, but involves no other unphysical quantities. In particular, the 
expression is independent of the box length L. 

The matrix element (L/7r)(0|r++(/i')|z) can be substituted directly 
to give an explicit expression for the two-point function. We see imme- 
diately that the correlator has the correct small-r behavior, for in that 
limit, it asymptotes to 




N^{2nb + nf) 



1 



which we expect for the theory of nh{nf) free bosons (fermions) at large 
K. On the other hand, the contribution to the correlator from strictly 
massless states is given by 



F(x ,x~^) 



E 



-(o|r++(i^)|i)" 

TT 



Mi=0 



It is important to notice that this l/r"^ behavior at large r is not the 
one we are looking for at large r. First of all, we do not expect any 
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massless physical bound state in this theory, and, additionally, it has 
the wrong N^. dependence. Relative to the 1/r^ behavior at small r, the 
behavior at large r that we expect is down by a factor of 1/Nc- This 
behavior is suppressed because we are performing a large- A^c calculation. 
All we can hope is to see the transition from the 1/r^ behavior at small 
r to the region where the correlator behaves like 1/r^. 

2.1.1 The A/" = (1, 1) theory. Although it is the J\f = (8,8) 
theory in which we arc ultimately interested in, wc can. nevertheless, 
perform the computation of the correlation function in models with less 
super symmetry. The evaluation of the correlator for the stress energy 
tensor in the M = (8, 8) theory is especially hard because of the many de- 
grees of freedom due to the large number of supercharges in that theory. 
We will show in Sec. 2.1.3 how to overcome this obstacle by exploiting 
a residual 'flavor' symmetry of the theory. To see how our numerical 
method works without these complications, it might be worthwhile to 
study the theory with supercharges (1,1). In the next section we will 
briefly cover the theory with a M = (2, 2) super symmetry. 

It has been argued that the J\f = (1,1) SYM theory does not exhibit 
dynamical supersymmetry breaking. A physicist's proof that supersym- 
metry is not spontaneously broken in this theory was given in Ref. [4]. 
This theory is also believed not to be confining [19] [20], and is therefore 
expected to exhibit non-trivial infra-red dynamics. The SDLCQ of the 
1+1 dimensional model with = (1, 1) supersymmetry was solved in 
Refs. [14, 21], and we apply these results directly in order to compute 
(23). For simplicity, we work at large Nc- The spectrum of this theory 
at finite K consists of 2K — 2 exactly massless states, i.e. K — I massless 
bosons, and their superpartners, accompanied by large numbers of mas- 
sive states separated by a gap. There is numerical evidence that this gap 
closes in the continuum limit. At finite we expect the degeneracy of 
2K — 2 exactly massless states to be broken, giving rise to precisely a 
continuum of states starting at M = as expected. 

The stress-energy correlator of this theory for various values of the 
harmonic resolution K, is shown in Fig. 3(a). We find the curious fea- 
ture that it asymptotes to the inverse power law c/r^ for large r. This 
behavior comes about due to the coupling (OlT^+Ira) with exactly mass- 
less states |n). The contribution to (23) from strictly massless states are 
given by 




L 



TT 




2 



87r2A;3 




Mn=0 



(24) 
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Figure 3. 



Log-Log plots of the two-dimensional correlation function 



{T++{x)T++{0)) (It) NH2u:tn,) ^-s- ^ at g^YMN/i: = L Left: (a) U = (1,1) 
SYM for A' = 4 to 10. Right: (b) JV = (2, 2) SYM for A' = 3, 4, 5, 6. 



TT 



-(0|r++(A:)|n) 



6 



M„=0 



(25) 



We have computed this quantity as a function of the inverse harmonic 
resolution 1/K and extrapolated to the continuum limit. The data cur- 
rently available suggests that the non-zero contribution from these mass- 
less states persists in this limit. 

2.1.2 The Af = (2, 2) theory. Let us now turn to the model 
with M = (2,2) supersymmetry. The SDLCQ version of this model 
was solved in Ref. [22]. The result of this computation can be inserted 
into Eq. (23). The result is shown in Fig. 3(b). This model appears to 
exhibit the onset of a gapless continuum of states more rapidly than the 
J\f = (1, 1) model as the harmonic resolution K is increased. Just as we 
found in the latter model, this theory contains exactly massless states 
in the spectrum. These massless states appear to couple to T^+jO) only 
for even K, and the overlap appears to be decreasing for growing K. It 
is believed that this model is likely to exhibit a power law behavior c/r^ 
for 7 > 4 for the T"^"^ correlator for r gYM\^ in the large Nc limit 
[4]. 

2.1.3 The Af = (8, 8) theory. In principle, we can now calcu- 
late the correlator numerically by evaluating Eq. (23). However, it turns 
out that even for very modest harmonic resolutions, we face a tremen- 
dous numerical task. At K = 2, 3, 4, the dimension of the associated 
Fock space is 256, 1632, and 29056, respectively. The usual procedure is 
to diagonalize the Hamiltonian P~ and then to evaluate the projection 
of each eigenfunction on the fundamental state r"'~+(— i^)|0). Since we 
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are only interested in states which have nonzero value of such projection, 
we are able to significantly reduce our numerical efforts. 

In the continuum limit, the result docs not depend on which of the 
eight supercharges Q~ one chooses. In DLCQ, however, the situation is 
a bit subtler: while the spectrum of (Q^)^ is the same for all a, the wave 
functions depend on the choice of supercharge [23]. This dependence is 
an artifact of the discretization and disappears in the continuum limit. 
What happens if we just pick one supercharge, say Since the state 
r++(— if )|0) is a singlet under R-symmetry acting on the "flavor" index 
of Q~, the correlator (23) does not depend on the choice of a even at 
finite resolution! 

We can exploit this fact to simplify our calculations. Consider an 

operator S commuting with both P~ and T'^~^{—K), and such that 
S\0) = so|0). Then the Hamiltonian and S can be diagonalized simulta- 
neously. We assume in the sequel that the set of states \i) is a result of 
such a diagonalization. In this case, only states satisfying the condition 
Sli) = sqIi) contribute to the sum in (23), and we only need to diagonal- 
ize P~ in this sector, which reduces the size of the problem immensely. 
We can deduce from the structure of the state T~^~^{—K)\0) that any 
transformation of the form 

ai^{k)^f{I)a^'^{k), /(/) = ±1 

6g(fe)^5(«)bg'"'(fc), ff(«)=±l (26) 

given arbitrary permutations P and Q of the 8 fiavor indices, commutes 
with T~^^{—K). The vacuum will then be an eigenstate of this trans- 
formation with eigenvalue 1. The requirement for P" = (Qi)'^ to be 
invariant under S imposes some restrictions on the permutations. In 
fact, we will require that Qi be invariant under S, in order to guarantee 
that P" is invariant. 

The form of the supercharge from [23] is 

Q-= [...]bUh)aj{ki)ai{k2) + ... (27) 
Jo 

+{l3i(5j - f3j(5j)ap[..]bl{k3)aj{ki)aj{k2) + .... 

Here the /?/ are 8x8 real matrices satisfying = 25jj. 

Let us consider the expression for Q^, Eq. (27). The first part of 
the supercharge does not include j3 matrices, and is therefore invariant 
under the transformation, Eq. (26), as long as g{l) = 1 and Q[l] = 1. We 
will consider only such transformations. The crucial observation for the 
analysis of the symmetries of the (3 terms is that in the representation of 

the P matrices we have chosen, the expression Bfj = (/3iPj— PjPJ) 
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may take only the values ±2 or zero. Besides, for any pair (/, J) there 
is only one (or no) value of a corresponding to nonzero B. Using this 
information, we may represent ;S in a compact form. With the definition 



together with the special choice of /? matrices we get the following 
expression for 
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The next step is to look for a subset of the transformations, Eq. (26), 
which satisfy the conditions g{l) = 1 and Q[l] = 1 and leave the matrix 
/X invariant. This invariance implies that 



The subset of transformations we are looking for forms a subgroup R 
of the permutation group Sg x Sg- Consequently, we will search for the 
elements of R that square to one. Products of such elements generate 
the whole group in the case of 5*8 x Ss- We will show later that this 
remains true for R. Not all of the Z2 symmetries satisfying (29) are 
independent. In particular, if a and b are two such symmetries then aba is 
also a valid Z2 symmetry. By going systematically through the different 
possibilities, we have found that there are 7 independent Z2 symmetries 
in the group R. They are listed in Table 1. We explicitly constructed 
all the symmetries of the type, Eq. (26), which satisfy Eq. (29) using 
Mathematica. It turns out that the group of such transformations has 
168 elements, and we have shown that all of them can be generated from 
the seven Z2 symmetries mentioned above. 

In our numerical algorithm we implemented the Z2 symmetries as 
follows. We can group the Fock states in classes and treat the whole 
class as a new state, because all states relevant for the correlator are 
singlets under the symmetry group R. As an example, consider the 
simplest non-trivial singlet 
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Table 1. Seven independent Z2 symmetries of the group R, which act on the 'flavor' 
quantum number of the difi'erent particles. Under the first of these symmetries, e.g., 
the boson oi is transformed into 07, etc. 



Hence, if we encounter the state a){l, l)a^(i^— 1, 1)|0) while constructing 
the basis, we will replace it by the class representative; in this case, by 
the state |1). Such a procedure significantly decreases the size of the 
basis, while keeping all the information necessary for calculating the 
correlator. In summary, this use of the discrete flavor symmetry of the 
problem reduces the size of the Fock space by orders of magnitude. 

In addition to these simplifications, one can further improve on the 
numerical efficiency by using Lanczos diagonalization techniques [25]. 
Namely, we substitute the explicit diagonalization with an efficient ap- 
proximation. The idea is to use a symmetry preserving (Lanczos) al- 
gorithm. If we start with a normalized vector \u\) proportional to the 
fundamental state r++(— if)|0), the Lanczos recursion will produce a 
tridiagonal representation of the Hamiltonian H^c = 2P'^P~ . Due to 
orthogonality of {|nj)}, only the (1,1) clement of the exponential of the 
tridiagonal matrix Hlc will contribute to the correlator [26]. We expo- 
nentiate by diagonalizing HlcVi = XiVi with eigenvalues Aj and get 

F(P^x^) = ^(J)'glte,,lV.&^ 

Finally, we Fourier transform to obtain 

which is equivalent to Eq. (23). This algorithm is correct only if the 
number of Lanczos iterations runs up to the rank of the original 
matrix, but in praxi already a basis of about 20 vectors covers all leading 
contributions to the correlator. 
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Numerical Results. To evaluate the expression for the correlator 
J'{r), we have to calculate the mass spectrum and insert it into Eq. (23). 
In the M = (8, 8) supersymmetric Yang-Mills theory the contribution of 
massless states becomes a problem. These states exist in the SDLCQ 
calculation, but are unphysical. It has be shown that theses states are 
not normalizable and that the number of partons in these states is even 
(odd) for K even (odd) [23]. Because the correlator is only sensitive 
to two particle contributions, the curves J^{r) are different for even and 
odd K. Unfortunately, the unphysical states yield also the typical 
behavior, but have a wrong Nc dependence. The regular 1/r^ contri- 
bution is down by 1/Nc, so we cannot see this contribution at large r, 
because we are working in the large Nc limit. 

We can use this information about the unphysical states, however, to 
determine when our approximation breaks down. It is the region where 
the unphysical massless states dominate the correlator sum. Unfortu- 
nately, this is also the region where we expect the true large-r behavior to 
dominate the correlator, if only the extra states were absent. In Fig. 4(a) 
for even resolution, the region where the correlator starts to behave like 
1/r^ at large r is clearly visible. In Fig. 4(b) we see that for even reso- 
lution the effect of the massless state on the derivative is felt at smaller 
values of r where the even resolution curves start to turn up. Another 
estimate of where this approximation breaks down, that gives consistent 
values, is the set of points where the even and odd resolution derivative 
curves cross. We do not expect these curves to cross on general grounds, 
based on work in [4], where we considered a number of other theories. 
Our calculation is consistent in the sense that this breakdown occurs at 
larger and larger r as K grows. 

We expect to approach the line dJ-{r)/dr = — 1 line signaling the 
cross-over from the trivial 1/r^ behavior to the characteristic 1/r^ be- 
havior of the supergravity correlator, Eq. (6). Indeed, the derivative 
curves in Fig. 4(b) arc approaching —1 as we increase the resolution and 
appear to be about 85 — 90% of this value before the approximation 
breaks down. There is, however, no indication of convergence yet; there- 
fore, we cannot claim a numerical proof of the Maldacena conjecture. A 
safe signature of equivalence of the field and string theories would be if 
the derivative curve would flatten out at —1 before the approximation 
breaks down. 

2.2. Three-dimensional correlators 

It remains a challenge to rigorously test the conjectured string/field 
theory correspondences. Although the so-called Maldacena conjecture 
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maybe the most exciting one, because it promises insight into full four- 
dimensional Yang- Mills theories in the strong coupling regime, there are 
other interesting scenarios. For instance, it was conjectured that the 
super gravity solutions corresponding to p + 1 SYM theories are black 
p-brane solutions, see e.g. Ref. [3]. Consequently, there are interesting 
testing scenarios also in three-dimensional spacetime. Numerically, of 
course, things get more difficult as the number of dimensions is increased. 
On the way to the full four-dimensional problem, it may be worthwhile 
to present our latest results on correlation functions in three dimensions; 
see also [13]. Fig. 7(b) shows the correlator for = 1 SYM(2+1) as 
a function of the distance r: it is converging well with the transverse 
cut-off T. To put things in perspective, we note that the construction 
of the largest Hamiltonian matrix involved in this calculation requires a 
Fock basis of approximately two million states. This is by a factor 100 
more than we used in the test of the Maldacena conjecture described in 
this article, which itself was already substantially better than the first 
feasibility study [4]. 

The correlator of the energy momentum operator has been studied in 
conformal field theory in 2-1-1 dimensions [27], and this provides a refer- 
ence point for our results. The structure of the correlators in conformal 
field theory is particularly simple in the collinear limit x± 0, and we 
therefore find it convenient to work in this limit. ^From results in con- 
formal field theory we expect that correlators behave as 1/r^ at small 
r, where we are probing deep inside the bound states. We have con- 
firmed this 1/r^ behavior by an analytic calculation of the free-particle 
correlator in the DLCQ formalism [16]. 

The contributions of individual bound states have a characteristic 
length scale corresponding to the size of the bound states. On dimen- 
sional grounds one can show that the power behavior of the correlators 
are reduced by one power of r; so for individual bound states the corre- 
lator behaves like 1/r^ for small r. It then becomes a nontrivial check 
to see that at small r the contributions of the bound states add up to 
give the expected l/r^ behavior. We find this expected result as well 
as the characteristic rapid convergence of SDLCQ at both small and 
intermediate values of r. 

At large r the correlator is controlled by the massless states of the 
theory. In this theory there are two types of massless states. At zero 
coupling all the states of the 1 + 1 dimensional theory are massless, and 
for non-vanishing coupling the massless states of the 1 + 1 theory are 
promoted to massless states of the 2 + 1 dimensional theory [17]. These 
states are BPS states and are exactly annihilated by one of the super- 
charges. This is perhaps the most interesting part of this calculation 
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because the BPS masses are protected by the exact supersymmetry of 
the numerical approximation and remain exactly zero at all couplings. 
Commonly in modern field theory one uses the BPS states to extrapo- 
late from weak coupling to strong coupling. While the masses of BPS 
states remain constant as functions of the coupling, their wave functions 
certainly do not. The calculation of the correlator at large r provides 
a window to the coupling dependence of the BPS wave functions. We 
find, however, that there is a critical coupling where the correlator goes 
to zero, which depends on the transverse resolution. A detailed study 
of this critical coupling shows that it goes to infinity linearly with the 
square root of the transverse resolution. Below the critical coupling the 
correlator converges rapidly at large r. One possible explanation is that 
this singular behavior signals the breakdown of the SDLCQ calculation 
for the BPS wave function at couplings larger than the critical coupling. 
If this is correct, calculation of the BPS wave function at stronger cou- 
plings requires higher transverse resolutions. We note that above the 
critical coupling (see Fig. 7 below) we do find convergence of the corre- 
lator at large r but at a significantly slower rate. 

Let us now return to the details of the calculation. We would like to 
compute a general expression of the form 

F(a;+,,T",,T^) = (0|r++(.x+,,T",.T^)r++(0,0,0)|0). (31) 

Here we will calculate the correlator in the collinear limit, that is, where 
X"*" = 0. We know from conformal field theory [27] calculations that this 
will produce a much simpler structure. 

The calculation is done by inserting a complete set of intermediate 
states |a), 

F{x+,x-,x^ = 0) = (32) 
^(0|r++(x-, 0, = 0)\a)e-'^^''^{a\T++{0, 0, 0)|0). 

a 

with energy eigenvalues P~. The momentum operator r++(x) is given 

by 



T++{x) = tr 



T++{x) + T++{x). 

(33) 

In terms of the mode operators, we find 



T++(x+,x-,0)|0) = -^^ ^ T(n,m)e-*(^"+^-)-nO), (34) 



2LI 

n,m n±,m± 
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-T++(n,m)|0) = ^tr 



and 



-ri+(n,m)|0) = — ^ — ^tr 

TT 4 



6t^(n,n_L)6]j(m,m_L) |0) . (36) 



Given each \a), the matrix elements in (32) can then be evaluated, and 
the sum computed. 

First, however, it is instructive to do the calculation where the states 
la) are a set of free particles with mass m. The boson contribution is 

2 



Fix+,x-,0)B= E ( 

n,m,s,t ^ 



4m 



(37) 



X V mnst (0|tr [a(n, n±)a{m, m±)]tT[a^ {s, s±)a^{t, t±)]\0), 
where the sum over n implies sums over both n and n±, and 



P- = 



m? + (2n_L7r//)' 
2n7r/L 



and P+ = ^ 



(38) 



The sums can be converted to integrals which can be explicitly evaluated, 
and we find 



F(x+,x-,0)b 



X 



2(27r)3 

where x"^ = 2x~x'^ . Similarly for the fermions we find 

/ vr ^2 



(39) 



F{x+,x-,Q)f= J2 

n,m,s,t 



\8m 



x(m - n)(s - t) (0|tr[6(n,n_L)fe(m,m_L)]tr[6l'(s,s_L)6l'(t,t_L)]|0). 
After doing the integrals we obtain 



F{x+,x-,0)f 



m" I X ' 



4(27r)3 X \x- 



K-r/2imx)Ks/2{mx) - Klj^imx) 



(41) 

We can continue to Euclidean space by taking r = \/2x+x to be real, 
and, finally, in the small-r limit we find 
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which exhibits the expected behavior. 

Now let us return to the calculation using the bound-state solution 
obtained from SDLCQ. It is convenient to write 

n,m,s,t ^ ' 

(43) 

where P^p is the Hamiltonian operator. We again insert a complete 
set of bound states \a) with light-cone energies P~ = {M'^ + P'j)/P^ 
at resolution K (and therefore P+ = ttK/L) and with total transverse 
momentum = 2N±'!r/l. We also define 

\u) = Nu-y2 ^n+m,KSnj_+m^,Nj_T{n,m)\0) , (44) 

TT — 

n,m 

where 7V„ is a normalization factor such that {u\u) = 1. It is straight- 
forward to calculate the normalization, and we find 

^ = ^(l-l)(2r + l). (45) 

The correlator (43) becomes 

2 



F{x+,x-,0)= (^) e-^''""^-^''^^"^|(^|a)r (46) 

We will calculate the matrix element {u\a) at fixed longitudinal res- 
olution K and transverse momentum N± = 0. Because of transverse 
boost invariancc the matrix element does not contain any explicit de- 
pendence on A^_|_ . To leading order in 1 /K the explicit dependence of the 
matrix element on K is K^; it also contains a factor of I, the transverse 
length scale. To separate these dependencies, we write F as 

F(x+x-Q)-— V 1 1 (^^y p-^+_^p+^- \{u\a)\^ 
^ ' ' ^ " 27r ^1^^^ 2LI\l) ' IKW ' ^ ^ 

Wc can now do the sums over K and A^_j_ as integrals over the longi- 
tudinal and transverse momentum components P^ = ttK/L and P"*" = 
27riVj_/Z. We obtain 



t) 



(48) 
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In practice, the full sum over a is approximated by a Lanczos [25] itera- 
tion technique [24, 26] that eliminates the need for full diagonalization of 
the Hamiltonian matrix. For the present case, the number of iterations 
required was on the order of 1000. 

Looking back at the calculation for the free particle, we see that there 
are two independent sums over transverse momentum, after the con- 
tractions are performed. One would expect that the transverse dimen- 
sion is controlled by the dimensional scale of the bound state {Rb) and 
therefore the correlation should scale like However, because of 

transverse boost invariance, the matrix element must be independent 
of the difference of the transverse momenta and therefore must scale as 



There are three commuting Z2 symmetries. One of them is the parity 
in the transverse direction, 

P : aij{k,n^) ^ -aij{k,-n^), bij{k,n^) ^ bij{k, -n^). (49) 

The second symmetry [28] is with respect to the operation 

S : aij{k,n'^) ^ —aji{k,n^), bij{k,n^) ^ —hji{k,n^). (50) 

Since P and S commute with each other, we need only one additional 
symmetry R = PS to close the group. Since Q~ , P and S commute 
with each other, we can diagonalize them simultaneously. This allows 
us to diagonalize the supercharge separately in the sectors with fixed P 
and S parities and thus reduce the size of matrices. Doing this one finds 
that the roles of P and S are different. While all the eigenvalues are 
usually broken into non-overlapping S-odd and S-even sectors [29], the 
P symmetry leads to a double degeneracy of massive states (in addition 
to the usual boson-fermion degeneracy due to supersymmetry). 

Numerical Results. The first important numerical test is the small- 
r behavior of the correlator. Physically we expect that at small r the 
bound states should behave as free particles, and therefore the correlator 
should have the behavior of the free particle correlator which goes like 
1/7.6 -^Ye see in (48) that the contributions of each of the bound states 
behaves like Therefore, to get the behavior of the free theory, 

the bound states must work in concert at small r. It is clear that this 
cannot work all the way down to r = in the numerical calculation. At 
very small r the most massive state allowed by the numerical approxi- 
mation will dominate, and the correlator must behave like To see 
what happens at slightly larger r it is useful to consider the behavior at 
small coupling. There, the larger masses go like 




(51) 
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Figure 5. The log-log plot of the correlation function / = 

r'5(T++(a;)r++(0)) (^f^) '^s- ^- ^eft: (a) in units where 

g = gYMVNJ/2Tr^/'^ = 0.10 for /s: = 4 and T = 1 to 9; Right: (b) in units 
where g = 5YM\/iVc7/27r^/^ = 1 for /sT = 5 and T = 1 to 9. 




Figure 6. Left: (a) Same as Fig. 5(b), but for g = gymV = 
10. Right: (b) Log-Log plot of the three-dimensional correlation function / = 

rHT++ix)T++{0)) (f^)' vs. r for p = gl.^Vm/27T'/' = 1.0 for = 6 

and T = 1 to 5. 

Consequently, as we remove the k'^ cutoff, i.e. increase the transverse 
resolution T, more and more massive bound states will contribute, and 
the dominant one will take over at smaller and smaller r leading to the 
expected This is exactly what we see happening in Fig. 5 at weak 

coupling with longitudinal resolution K = 4 and 5. 

The correlator converges from below at small r with increasing T, 
and in the region —0.5 < logr < 0.5 the plot of times the correlator 
falls like 1/r. In Fig. 7 at resolution = 5 we see the same behav- 
ior for strong coupling {g = 5YM\/A^/27r3/2 = 10) but now at smaller r 
(log r c± —0.5) as one would expect. Again at strong coupling we see that 
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the correlator converges quickly and from below in T. All indications 
are that at small r the correlators are well approximated by SDLCQ, 
converge rapidly, and show the behavior that one would expect on gen- 
eral physical grounds. This gives us confidence to go on to investigate 
the behavior at large r. 

The behavior for large r is governed by the massless states. From 
earlier work [30, 17] on the spectrum of this theory we know that there 
are two types of massless states. At g' = the massless states arc a 
reflection of all the states of the dimensionally reduced theory in 1 + 1. 
In 2+1 dimensions these states behave as g'^M-^^^. We expect therefore 
that for g — there should be no dependence of the correlator on the 
transverse momentum cutoff T at large r. In Fig. 5(a) this behavior is 
clearly evident. 

At all couplings there are exactly massless states which are the BPS 
states of this theory, which has zero central charge. These states are 
destroyed by one supercharge, Q~, and not the other, Q"*". From ear- 
lier work [30] on the spectrum we saw that the number of BPS states 
is independent of the transverse resolution and equal to 2K — 1. Since 
these states are exactly massless at all resolutions, transverse and lon- 
gitudinal convergence of these states cannot be investigated using the 
spectrum. These states do have a complicated dependence on the cou- 
pling g through their wave function, however. This is a feature so far 
not encountered in DLCQ [16]. In previous DLCQ calculations one al- 
ways looked to the convergence of the spectrum as a measure of the 
convergence of the numerical calculation. Here we see that it is the cor- 
relator at large r that provides a window to study the convergence of the 
wave functions of the BPS states. In Fig. 7 we see that the correlator 
converges from above at large r as we increase T. 

We also note that the correlator at large r is significantly smaller than 
at small r, particularly at strong coupling. In our initial study of the 
BPS states [17] we found that at strong coupling the average number 
of particles in these BPS states is large. Therefore the two particle 
components, which are the only components the T~^~^ correlator sees, 
are small. 

The coupling dependence of the largc-r limit of the correlator is much 
more interesting than we would have expected based on our previous 
work on the spectrum. To see this behavior we study the large-r behavior 
of the correlator at fixed g as a function of the transverse resolution T 
and at fixed T as a function of the coupling g. We see a hint that 
something unusual is occurring in Fig. 7. For values of the coupling up 
to about g = 1 we see the typical rapid convergence in the transverse 
momentum cutoff; however, at larger coupling the convergence appears 
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Figure 7. The large-r limit of the log of the correlation function / = 
r^T++{x)T++{0)) (^^y if^^^ vs. 1/T for [left](a) X = 5 and [right](b) K = 6 
and for various values of the coupling g = gYMVNcl/2n^^'^. 




Figure 8. The largo-r limit of the correlation function / = 

r^(T++(x)T++(0))(f^)'^^ vs. g = gYMVKi/2^'/' for [left](a) K = 5 
[right] (b) K = 6 and for various values of the transverse resolution T. 



to deteriorate, and we see that for g = 5 the correlator is smaller than 
Ski g = 10. We see this same behavior at both K = 5 and K = 6. We do 
not see this behavior at K = A, but it is not unusual for effects to appear 
only at a large enough resolution in SDLCQ. In Fig. 8 we see that the 
correlator does not in fact decrease monotonically with g but rather has 
a singularity at a particular value of the coupling which is a function of 
K and T. Beyond the singularity the correlator again appears to behave 
well. 

If we plot the 'critical' couplings, at which the correlator goes to 
zero, versus VT, as in Fig. 9, we see that they lie on a straight line, 
i.e. this coupling is a linear function of s/T in both cases, K = 5 and 
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Figure 9. Critical coupling gcrit versus VT for (a) K = 5 and (b) K = 6. 



6. Consequently, the 'critical' coupling goes to infinity in the transverse 
continuum limit. It appears as though we have encountered a finite 
transverse cutoff effect. The most likely conclusion is that our numerical 
calculation of the BPS wave function is only valid for g < gcnt{T). 
While the large-r correlator does converge above the critical coupling, 
it is unclear at this time if it has any significance. It might have been 
expected that one would need larger and larger transverse resolution to 
probe the strong coupling region, the occurrence of the singular behavior 
that we see is a surprise, and we have no detailed explanation for it 
at this time. We see no evidence of a singular behavior at small or 
intermediate r. This indicates, but does not prove, that our calculations 
of the massive bound states is valid at all g. 

We do not seem to see a region dominated by the massive bound 
states, that is, a region where r is large enough that we see the structure 
of the bound states but small enough that the correlator is not dominated 
by the massless states of the theory. Such a region might give us other 
important information about this theory. 

3. Conclusions 

In this note we reported on progress in an attempt to rigorously test 
the conjectured equivalence of two-dimensional = (8, 8) supersym- 
metric Yang- Mills theory and a system of Z)l-branes in string theory. 
Within a well-defined non-perturbative calculation, we obtained results 
that are within 10-15% of results expected from the Maldacena conjec- 
ture. The results are still not conclusive, but they definitely point in 
right direction. Compared to previous work [4], we included orders of 
magnitude more states in our calculation and thus greatly improved the 
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testing conditions. We remark that improvements of the code and the 
numerical method are possible and under way. During the calculation 
we noticed that contributions to the correlator come from only a small 
number of terms. An analytic understanding of this phenomenon would 
greatly accelerate calculations. We point out that in principle we could 
study the proper 1/r behavior at large r by computing 1/Nc corrections, 
but this interesting calculation would mean a huge numerical effort. 

In this work we also discussed the calculation of the stress-tensor 
correlator (0|r++(x)r++(0)|0) in A/" = 1 SYM in 2+1 dimension at large 
Nc in the collincar limit. We find that the free-particle correlator behaves 
like in agreement with results from conformal field theory. The 

contribution from an individual bound state is found to behave like 
and at small r such contributions conspire to reproduce the conformal 
field theory result We do not seem to find an intermediate region 

in r where the correlator behaves as refiecting the behavior of the 
individual massive bound states. 

At large r the correlator is dominated by the massless BPS states of 
the theory. We find that as a function of g the large-r correlator has a 
critical value of g where it abruptly drops to zero. We have investigated 
this singular behavior and find that at fixed longitudinal resolution the 
critical coupling grows linearly with \/T. We conjecture that this critical 
coupling signals the breakdown of SDLCQ at sufficiently strong coupling 
at fixed transverse resolution, T. While this might not be surprising 
in general, it is surprising that the behavior appears in the BPS wave 
functions and that we see no sign of this behavior in the massive states. 
We find that above the critical coupling the correlator still converges but 
significantly slower. It is unclear at this time if we should attach any 
significance to the correlator in this region. 

This calculation emphasizes the importance of BPS wave functions 
which carry important coupling dependence, even though the mass eigen- 
values are independent of the coupling. We will discuss the spectrum, 
the wavefunctions and associated properties of all the low energy bound 
states of A/" = 1 SYM in 2-1-1 dimensions in a subsequent paper [31]. 

A number of computational improvements have been implemented in 
our code to allow us to make these detailed calculations. The code now 
fully utilizes the three known discrete symmetries of the theory, namely 
supersymmetry, transverse parity P, Eq. (49), and the Z2 symmetry 
S", Eq. (50). This reduces the dimension of the Hamiltonian matrix 
by a factor of 8. Other, more efficient storage techniques allow us to 
handle on the order of 2,000,000 states in this calculation, which has 
been performed on a single processor Linux workstation. Our improved 
storage techniques should allow us to expand this calculation to include 
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higher supersymmetries without a significant expansion of the code or 
computational power. We remain hopeful that porting to a parallel 
machine will allow us to tackle problems in full 3+1 dimensions. 
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